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Abstract. We apply a novel technique of Faraday Rotation measure (RM) map analysis to three galaxy clusters, Abell 400, 
Abell 2634 and Hydra A, in order to estimate cluster magnetic field strengths, length scales and power spectra. This analysis - 
essentially a correlation analysis - is based on the assumption that the magnetic fields are statistically isotropically distributed 
across the Faraday screen. We investigate the difficulties involved in the application of the analysis to observational data. We 
derive magnetic power spectra for three clusters, i.e. Abell 400, Abell 2634 and Hydra A, and discuss influences on their shapes 
caused by the observational nature of the data such as limited source size and resolution. We successfully apply various tests 
to validate our assumptions. We show that magnetic fluctuations are probed on length scales ranging over at least one order 
of magnitude. Using this range for the determination of magnetic field strength of the central cluster gas yields 3 /jG in Abell 
2634, 6 yuG in Abell 400 and 12 /jG in Hydra A as conservative estimates. The magnetic field autocorrelation length Ab was 
determined to be 4.9 kpc for Abell 2634, 3.6 kpc for Abell 400 and 0.9 kpc for Hydra A. We show that the RM autocorrelation 
length /{rm is larger than the magnetic field autocorrelation length Ag - for the three clusters studied, we found /Irm ^ L.AAg 
- and thus, they are not equal as often assumed in the literature. Furthermore, we investigate in a response analysis if it is 
possible to determine spectral slopes of the power spectra. We find that integrated numbers can be reliably determined from 
this analysis but differential parameters such as spectral slopes have to be treated differently. However, our response analysis 
results in spectral slopes of the power spectra of spectral indices a = 1.6 to 2.0 suggesting that Kolmogorov spectra (a = 5/3) 
are possible but flatter spectral slopes than o- = 1.3 can be excluded. 

Key words. Radiation mechanism: non-thermal - Galaxies: active ~ Intergalactic medium - Galaxies: cluster: general - Radio 
continuum: general 



1 . Introduction iBmnetti et alJbOoH) . Another possibility is the comparison of 

their synchrotron radiation fluxes with inverse Compton X- 

^ Clusters of galaxies are inhabited by magnetised plasma. The ^^^^^on leading to lower limits to the field strengths of 

detection of diff^use synchrotron emission observed as cluster ^^^^^ 0.2...1 //G (e.g.jar ris & Grindl a^fl97ilRe phaeli 

radio halos and radio relics is but one piece of evidence for r^. ^ Bie rmannlll998t iFusco-Femiano et alJll9^ 

the existence of clu s ter wide magnetic fields (for a review, ,Bagchi et al.l99&) 
see ICarilH & Tavio3200i) . Still, their strengths, their structures 

and their origin remains to be unveiled. A different method to gain knowledge about magnetic fields 

Although the cluster wide magnetic fields are not believed the observation of the Faraday rotation effect which arises 

to influence the dynamics directly in most clusters because flie whenever polarised radio emission passes through a magne- 

magnetic pressure is estimated to be too small compared to the ^ised medium, causing its plane of polai-isation to rotate if the 

thermal pressure, they are still important for the understanding magnetic field component parallel to flie direction of the light 

of die physics of flie intra-cluster gas. Indirectly fliey ai-e likely propagation is nonzero. If the Faraday active medium is exter- 

an important cluster gas ingredient, e.g. cluster magnetic fields ^o the source of the polarised radio emission one expects 

suppress heat conduction as observed in clust er cold fronts (e.g. *e rotation to be proportional to flie squared wavelengfli. The 

iMarkevitch et alJ200(Alvikhiinin et alJ200ll) . proportionality constant is called the Faraday rotation measure 

rru A-cc . »i, J » ■ ■ c u . (RM) which can be evaluated in terms of the line of sight inte- 

There are different methods to gam information about ^ ' ° 

^. c , . ^ . ^ ^, , ^ .. ^. gral over the product of the electron density and the magnetic 

magnetic fields. One is to use the synchrotron radiation ^ t- j b 

, J J- u 1 1- » » *u ij field component along the fine of sight, 

observed as radio halos or relics to estimate the field ^ ^ ° 

strength. Energy equipartition arguments for these clusters Recently, 'Rudnick & Blundelf ('2003') claimed to have ev- 

lead to strengths of about ~ 0.1. ..2 yG (e.g. iFerettL .1999^ idence for the RM producing fields being local to the po- 
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larised radio source as suggested by iBicknell et al.' (!l99d). 
However, this conclusion was based on a biased statistic with- 
out performing a proper null-experiment (Enfi lin et al.l l2003). 
Furthermore, the assumption that the Faraday rotation is pro- 
duced in the intra-cluster gas is also supported by indepen- 
dent evidence. We mentioned already some of this evidence, 
the cluster wide diffuse radio emission observed as radio halos 
and relics. Further evidence is the observation of the Laing- 
Garrington effect. It manifests itself by an asymmetric Faraday 
depolarisation of double radio lobes located within galaxy clus- 
ters (Garrington et al. 1988; Laing 1988). This asymmetry is 
caused by different path lengths for radio emission through the 
Faraday active m edium between source and observer as sug- 
gested by Garring ton & Conwavl ( ll99ll) . This effect could also 
be explained by an asymmetric mixing layer being thin close to 
the head side of a FR II while getting thicker on the back-flow 
side of it. However, this explanation fails to give a convincing 
explanation for the large RM and dep olarisation asymmetry ob- 
served for the FR I source Hydra A ( Tavlor & Perlev 1993) . 

Further independent evidence for cluster wide magnetic 
fields is provided by a rece nt statistical RM inv estigation of 
point sources performed by IClarke et al.l (l200ll) . It revealed 
larger RM values for sources observed through the intra-cluster 
gas in comparison to a control sample of sources where no 
intra-cluster gas was located between source and observer This 
suggests that the enhancement of the RM towards the cluster 
centre results most probably from the magnetised cluster gas. 

Thus, we believe that the analysis of Faraday rotation maps 
of extended polarised radio sources located behind or embed- 
ded in galaxy clusters allows us to understand the strength and 
to get some hints about the structure of cluster magnetic fields. 

A naive analysis of such maps is to estimate a minimum 
field strength by assuming a constant magnetic field through- 
out the cluster along the line of sight. However, close inspection 
of RM maps of extended extragalactic radio sources reveals a 
patchy structure on small scales (5 - 20 kpc) which indicates 
that the fields are tangled. An improved analysis of the RM 
maps assumes the Faraday screen consists of cells which have 
a constant size and a constant magnetic field strength but the 
field direction varies randomly from cell to cell. The mean RM 
which would be produced by such a screen builds up in a ran- 
dom walk and thus has a zero mean but a non vanishing disper- 
sion. This dispersion will be proportional to the square root of 
the number of cells along the fine of sight, the cell size, the elec- 
tro n densi ty profile a nd the magnetic fie ld strength as shown 
bv lFeltenI (f 19961) and lFeretti et all ( Il995h . Using these propor- 
tionalities magnetic field strengths of the order of a few iiG 
for no n-cooling flo\y clus ters, e.g. Coma ( Feretti et al. 1995), 
A119 (iFe retti et al Jll999l) . 3C129 llTavlor et alJl2001l) A514 
mietal...20Qll) . A400 & A2634 feilek & Owen"200^, 
and of the order of a few 1 /iG for cooling flow clusters, e.g. 
Cygnus A (Dreheret al. 1987), 3C295 (Perley & Taylor 1991), 
A 1795 (.Ge & Owea.1993.) . Hydra A OTavlor & Perley..l993) . 
were derived. 

However, there are some drawbacks of the cell model. First 
of all, the divergence of the magnetic field in such a Faraday 
screen would be nonzero which contradicts Maxwell's equa- 



tions. Secondly, a spectrum of scales for the cell sizes is more 
likely than a single scale. 

Another important issue in previous analyses is the assump- 
tion that the RM ordering scale read off RM maps is equivalent 
to the magnetic field's characteristic length scale, the autocor- 
relation length Ab- This might have lead to underestimations of 
field strengths in the past since the equation often used for the 
derivation of the magnetic field strength is 



(RM^) 
nl L A- 



(1) 



where (RM ) is the RM dispersion, L is the depth of the 
Faraday screen, aq - l{2nm^gC^), is the electron density 
and /i, is a characteristic length scale of the fields. At this point 
it becomes clear that the definition of the characteristic length 
scale /i, is crucial for the derivation of magnetic field strength. 
The assumption of a constant magnetic field throughout the 
cluster leads to a definition of A- - L and thus, (B?) oc L'^. 
The characteristic length scale for the cell model would be the 
size of each cell A- - Iceii- Another definition one could think 
of is the RM ordering scale. A- - /Irm- However, the correct 
length scale is the autocorrelation length A~ of the z-component 
of the magnetic field measured along the line of sight. 

A purely statistical approach which incorporates the van- 
is hing divergence of the magnetic field was recently developed 
bv lEnBlin & Vo gt (2003). This approach relies on the assump- 
tion that the fields are statistically isotropically distributed in 
Faraday screens such that the z-component is representative for 
all components. Starting from this assumption, a relation be- 
tween the observationally accessible RM autocorrelation func- 
tion and the magnetic autocorrelation tensor is established in 
real and in Fourier space such that one gains access to the 
power spectrum of the magnetic field inhabited by the Faraday 
scr een. 

lEnBlin & Vog tl (|2003 ) demonstrated that in the case of 
isotropy the magnetic autocorrelation length A^ in Eq.Q is 
given by A, = 3/2 Ab, where Ab is the 3-dimensional magnetic 
autocorrelation length. Therefore, Eq. Q becomes 



2<RM2) 



LA, 



2 7 
0(1 



(2) 



The length scale Ab can be estimated from the measured RM 
power spectrum. 

Here we apply this statistical analysis to observational data 
by reanalysing the Faraday rotation measure maps of three ex- 
tended extragalactic radio sources: Hydra A ('Ta ylor & Perlevl 
1993), 3C75 and 3C465 iEilek & Owen 2002) which were 
kindly provided by Greg Taylor, Frazer Owen, and Jean Eilek, 
respectively. 

In section |2l we outline shortly the idea of our approach, 
the real and the Fourier space analysis and what possibili- 
ties are given to test if the assumptions made are reasonable. 
Furthermore, we state the equations which were used in the 
analysis wh ile their derivation and more detailed discussion can 
be found in EnBlin & Vogt (2003). In section|31 the application 
of these approaches to the observational data and the difficulties 
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involved are described. Emphasis is given to a critical discus- 
sion of the strengths and limits of the approach. The results are 
presented and discussed in section 0] Finally, the conclusions 
are drawn in the last section|5] 

Throughout the rest of the paper we assume a Hubble con- 
stant of Ho = 70kms"' Mpc ',Q,„ = 0.3 and Q a ^ 0.7 in a flat 
univer se. All equations follow the notation of lEnBlin & Vog tl 
ll2003l) . 

2. Method 

2.1. Idea 

We assume that the magnetic fields in galaxy clusters are 
isotropically distributed throughout the Faraday screen. If one 
samples these field distributions over a sufficiently large vol- 
ume they can be treated as statistically homogeneous and sta- 
tistically isotropic. This means that the statistical average over 
any field quantity will not be influenced by the geometry or the 
exact location of the volume sampled. 

One of the magnetic field quantities which is of inter- 
est here is its autocorrelation function since this function can 
be used to calculate its correlation length and field strength. 
Furthermore, the autocorrelation function is equivalent to the 
Fourier transform of the power spectrum as stated by the 
Wiener- Khinchin Theorem (WKT). 

Since the magnetic field is a vector quantity, an autocorre- 
lation tensor has to be considered rather than a function. The 
diagonal elements of this tensor are all equal in the case of 
magnetic isotropy which we assume as stated above. On the 
other hand, its off-diagonal terms consist of two numbers de- 
scribing the symmetric and non-symmetric contribution. The 
divergence-freeness of magnetic fields (V • B = 0) couples 
the symmetric components of the tensor which means that the 
knowledge of a diagonal element ensures complete knowledge 
of the symmetric part of the tensor Introducing now the scalar 
magnetic autocorrelation function w{r) as the trace of the au- 
tocorrelation tensor w{r) - {B(x) ■ Six + r)) i, field quantities 
can be derived such as the magnetic autocorrelation length and 
average magnetic field energies using this function. 

The non-symmetric part of the tensor represents the mag- 
netic helicity. These off'-diagonal terms of the tensor do not en- 
ter in the definition of the scalar magnetic autocorrelation func- 
tion w{r) and thus, do not have any effect on the quantities of 
interest listed above. 

The observable considered in our approach is the Faraday 
rotation which is related to the projected magnetic field distri- 
bution along the line of sight. It is possible to establish a rela- 
tion between the autocorrelation function of the RM map and 
the autocorrelation function of the magnetic field as discussed 
in more detail in the next sections. 

As stated above, the knowledge of the autocorrelation 
function gives access to the power spectra using the WKT. 
However, any phase information is lost by measuring power 
spectra. Therefore, it is possible that different realisations of 
magnetic field structure have the same power spectrum. This 
implies that not all information which might be contained in the 
RM maps are accessible by the analysis presented. Correlation 
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function of higher order statistic are needed to access these in- 
formation. 

In the discussion so far, we ignored the fact that we are us- 
ing observational data and that the RM is also dependent on 
electron density and magnetic energy density on global scales 
which makes the introduction of a window function necessary. 
This window function can be interpreted as the sampling vol- 
ume. It accounts for the limited size of the radio source and the 
global properties of the cluster gas such as the electron density 
and the average magnetic energy density profiles. A virtually 
statistical homogeneous magnetic field can be thought to be 
observed through this function. 

The window function will be zero for points outside the 
volume which is situated in front of the radio source. It will 
scale with the electron density and the global magnetic field 
energy distribution within the Faraday screen between radio 
source and observer. Furthermore if wanted the window func- 
tion can contain noise reducing data weighting schemes. 

The overall effect of any finite window in real space is to 
smear out power in Fourier space. Furthermore, any analysis 
employing a window function will be sensitive to magnetic en- 
ergies on scales smaller than the window size and insensitive 
to scales larger than the window size. Therefore, understanding 
the influence of the window function is crucial for the evalua- 
tion of the results obtained by applying of our analysis to ob- 
servational data. Possibilities to assess its influence are given 
in section lT^ 

2.2. Real Space Analysis 

For a line of sight parallel to the z-axis and displaced by x± 
from it, the Faraday rotation measure arising from polarised 
emission passing from the source at Zfi(x_i_) through a magne- 
tised medium to the observer located at infinity is expressed 
by 

/-»oo 

RMixJ^aol dznAx)B,(x), (3) 

where - e-'/(27rmgC^), x = ix±,z), ng(x) is the electron 
density and B.{x) is the magnetic field component parallel to 
the line of sight. 

The autocorrelation function of the Faraday rotation map 
can be defined as 

CRM(r^) = <RM(;«:^)RM(;«:^ + rj},^ , (4) 

where the brackets mean the map average with respect to x±. 
However, an observed RM map is limited in size which leads 
to noise in the calculation of the correlation function on large 
scales since less pixel pairs contribute to the correlation func- 
tion for larger pixel separations. For suppressing this under- 
sampling, the observable correlation function was chosen to be 

CrA(>-±) = ^ / ^-^i RM(^±)RM(;t:^ + rj, (5) 

where the area An of the region Q for which RM's are actually 
measured is used as normalisation assuming that RM(jrj.) = 
for x± i Q.. Furthermore, it is required that the mean RM is 
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zero since statistically isotropic divergence free fields have a 
mean RM of zero to high accuracy. The non vanishing mean 
RM in the observational data stems from foregrounds (e.g. the 
galaxy) in which we are not interested or has its origin in large 
scale fields to which the approach is insensitive by construction 
in order to suppress statistically under-sampled length scales. 

As stated in section ITTI the virtually homogeneous mag- 
netic field component can be thought to be observed through 
a window function f{x) which describes the sampling volume. 
One would chose a typical position in the cluster jCief (e.g. its 
centre) and define n^o - ndXref) and Bo - (^^(J^ref))'''^- The 
window function can then be expressed as 



fix) = l|;c^en| 1|z>z,(a:^)| h(xj_) g(x) n,(x) /n,o. 



(6) 



where l{condition) is equal to one if the condition is true and 
zero if not. The dimensionless average magnetic field profile 
g(x) = (B^(jc))'^^/Bo is assumed to scale with the density 
profile such that g(x) - (iieix) I Heof" . Reasonable values for 
the exponent an r ange between 0.5 and 1 as suggested by 
iDolag et'aP ( 1200 ih . The function h{x) allows to introduce a 
suitable data weighting scheme and can be chosen to be unity 
if no data weighting is required. 

As mentioned above, the magnetic autocorrelation function 
H'(r) can be defined as the trace of the magnetic field auto- 
correlation tensor M,j = {Biix) Bj(x + r )) assu ming isotropy 
w(r) = {B(x) ■ B(x + r)} = Zi Ma, re.g. lSubramanian 1999 ). 
However, the Faraday rotation is also related to the magnetic 
field via Eq.Q- Therefore any RM correlation contains also 
information on the autocorrelation of the magnetic field com- 
ponent along the line of sight 



dr.M,,Sr) 

oo 

^ J — OO 



(7) 
(8) 



where r - {r^,rS). Taking the effect of the window function 
into account, we introduced C±{r±) - {C'^^(r±))/au where 
ai = OqIi^qL. The parameter L can be interpreted as the char- 
acteristic depth of the Faraday screen and can be expressed 
as L = V[f]/AQ, where V[f] is the probed efifective volume 
Jdx'fix). 

Using Eq. one is now able to establish a relationship 
between the RM autocorrelation function C±(r±) and the mag- 
netic autocorrelation function w(r). This relationship can be ex- 
pressed as 



w(r) 



2 d r 

nr dr J,. 
= -t dy 



' yCAy) 

dy — 

C'(y) 

^y2 _ ^2 



(9) 



(10) 



where for the deprojection of Eq. an Abel integral equation 
was used assuming w(r) stays bounded for r — > oo. 

These two equations enable us to calculate not only the 
magnetic autocorrelation function but also to obtain directly 
the average magnetic energy density (sb) by using the relation 
(bb) - w(0)/87r and thus, we are able to determine the magnetic 



field strength by calculating the value of the magnetic autocor- 
relation function at point r = 0: w{r = 0) = {B^}. 

The RM autocorrelation length /Irm and the magnetic auto- 
correlation length Ab can be calculated by integrating the cor- 
relation functions: 



•M - I 



dr 



dr_i_ 



<RM2(a:)) 

XOO 
_^drrw{r) 



dr 



CAQ) r drw(r) 

J— CO ^ 

{B(x + re)-B(x)h,e 
(BHx)) 



, w(r) C^(0) 

dr = 2 , 

w(0) w(0) 



(11) 
(12) 

(13) 
(14) 



where e and Cj. are 3- and 2-dimensional unit vectors, respec- 
tively, over which we are averaging. 

2.3. Fourier Space Analysis 

In Fourier space a similar formalism can be developed assum- 
ing statistical isotropy and relying on the V • B - condition. 
After performing a 2-dim Fourier transform' of the RM map, 
one can calculate the RM autocorrelation in Fourier space: 



CJkJ 



<|RM(fcJp) 
aiAn 



(15) 



implied by the WKT which connects the 2-dimensional Fourier 
transformed RM maps with the RM autocoiTelation function. 

In Fourier space Eq.(|8}, the relationship between the RM 
and the magnetic field autocorrelation function, can be ex- 
pressed by 

CAk±) = Mu.(k^,0) = ^w(k^,0). (16) 

This equation states that the k^- = plane of the M-r component 
(see Eq. (0) of the Fourier transformed magnetic autocoiTela- 
tion tensor is accessible by the 2-dimensional Fourier transfor- 
mation of the RM map. Since all diagonal elements of the au- 
tocorrelation tensor are equal in the case of statistical isotropy, 
no more information is required in order to reconstruct the full 
magnetic autocorrelation. 

The 3-dimensional power spectrum w{k)d^k is related to 
the 1 -dimensional magnetic energy spectrum in the Fourier 
space sf'(k) dk by 



2(27r) 



3 ' 



(17) 



where w(k) = w(k) because of isotropy. 

The combination of Eq. ( I15t and Eq. leads to an equa- 
tion which relates the magnetic energy spectrum directly to the 
Fourier transformed RM map 



sf^k) 



In 



fl'0|RM(^^)|2 



(18) 



' defining the Fourier transformation for a n-dimensional function 
F{x) such that: 

F{k) = / d"xFix)e''' '';F{x) = f d"kF{k)e-''' ''. 
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which is equivalent to averaging over rings in k-space. 
Therefore integration over k yields the average magnetic en- 
ergy density s°g^ - dks°^^(k) and thus, leads to a magnetic 
field strength estimate. 

Other quantities like the correlation length scales, Ag and 
Arm, can also be calculated via integrating over the correlation 
function: 



= 2 



rdkw(k) 
Jg dk k w(k) 
j^dkkwik) 
^ J^^dkk^ w(k)' 



(19) 
(20) 



From these two equations, it is obvious that the RM auto- 
correlation length /Irm has a much larger weight on large-scale 
fluctuations than the magnetic correlation length Ag- Therefore 
one expects that /Irm > Ab in the natural case of a broadly 
peaked magnetic power spectrum. 

The autocorrelation function and the power spectrum are 
related by a Fourier transform as stated by the WKT, therefore 
one can derive a relation for the magnetic autocorrelation func- 
tion w{r) in terms of the power spectrum w{k) 



w(r) 



" W Jo 



dkk^ w(ky 



sin kr 



kr 



(21) 



The derivation of the function w(r) employing the power spec- 
tra wiK) can be used as a good test for proper numerics, since 
it can be compared to the deprojected magnetic autocorrelation 
function w{r) determined in real space. 

2.4. Testing the model and assessing tlie influence of 
the window function 

There are dififerent methods to assess the results obtained us- 
ing the approaches in real and Fourier space described above. 
Besides the comparison of the numbers determined for the 
characteristic length scales and magnetic field strengths in both 
spaces, there is also the possibility to evaluate the quality of 
the models used for the magnetic energy density and the loca- 
tion of the radio emitting source within or behind the Faraday 
screen. For this purpose a ;t'^-test function can be introduced 
such that 



2 RM(xJ' 



(22) 



where the expected RM dispersion is defined as {RM(acj^)^) = 
0.5 aQn^QEQAs J_^dzf^(x). If one averages this function radi- 
ally and examines it for any tendency in its spatial distribu- 
tion, one has a first test of the model. If there are any large 
scale trends observed one should reconsider the model used. 
Furthermore, the integration of this function should result in a 
value close to unity 



An I 



dx]_x^(xj_) ~ 1, 



(23) 



if the choice of the model is a suitable description of the global 
magnetic energy distribution and the geometry of the source. 



Higher values of xly would arise for unsuited models and one 
would have to refine the model parameters. 

The isotropy which is an important assumption for our anal- 
ysis can be assessed by the inspection of the Fourier trans- 
formed RM maps. The assumption of isotropy is valid in the 
case of a spherically symmetric distribution of |RM(A;)p. 

A completely different problem is the assessment of the in- 
fluence of the window function on the shape of the magnetic 
power spectrum and thus, on the measured magnetic energy 
spectrum. As stated above, the introduction of a finite window 
function has the effect to smear out the power spectrum and 
Ea.( ll6t becomes 



d q w{q) -J 



ql \f(ki. - qr 



(24) 



The term q]_/q^ < 1 in the expression above states some loss 
of magnetic power. Without this term, the expression would 
describe the redistribution of magnetic power within Fourier 
space, where in this case the magnetic energy would be con- 
served. 

However, the expression can be employed as an estimator 
of the response of the observation to the magnetic power on a 
given scale p by setting w{q) - 6{q - p) in the above equa- 
tion. For an approximate treatment of a realistic window and a 
spherical data average, one can derive 



2p r>"^'?i; 

= — \dq 



d(pW±(q(cos (p, sin <;*)) 



(25) 



{q^ + p^ 



■klf 



(26) 



Here the projected Fourier window was introduced 

(27r)2A[/j 

where A[/j^] = J t/^x^ /^(jCj_) while defining the projected win- 
dow function as /^(jri) = ^ dzf^{x)VYf]IAn- 

One can also require for the input power spectrum a field 
strength by setting w{q) - In^ (B^ / p^) 6{q - p) in Eq. 
The resulting power spectrum can be treated as an observed 
power spectrum and thus, a magnetic field strength Z?exp can be 
derived by integration following Eq. (I17> . The comparison of B 
and Bexp for the different p-scales yields a further estimate for 
reliable ranges in k-space. 

So far the response power spectrum was derived for an 
underlying magnetic spectrum located only on scales p. One 
could also add the particular response functions such that the 
actual observed power spectrum w°'°^{k) would be matched. 
Typical power spectra, e.g. for turbulence, show a power 
law behaviour over several orders of magnitude. Therefore, 
one would choose to weight the different response functions 
w^p'^{ki_) derived using Eq. (125 > by a power law while integrat- 
ing over them 



-a-2 



(27) 



where wq represents a normalisation constant, fco a typical k- 
space length scale and pmin and pmax are upper and lower cut- 
offs for the integration over the response functions. Ideally, one 
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would choose the normalisation constant wq = cq B^^jk^, where 
Co is chosen such that ^ dk^ w'^'^'^{k) - Bq. One would vary the 
spectral index a of the power law, Bo or wo and the lower cut- 
off of the integration in order to match the two functions, 
w°'''(k) and w'^^'^ik). 

One has direct access to the average magnetic energy den- 
sity sb by integrating Ea.j27ll which results in the analytic ex- 
pression 



SB 



Wo kl 



2(27r)3(Q'- 1) 



^0 



(28) 



Following the definition of the magnetic autocorrelation 
length Ab (Eq. i2Q\ ). one can derive an expression for this 
length using the parameters determined in the analysis above 



1 ^'min Pm-dx 



^mm /-"max 



(29) 



In the calculation of the expected response C^^^{k^) and 
thus, M>™P(^x) to an underlying magnetic power spectrum w(k) 



seen through some win dow W(k) as exp ressed by Eq. (I25> (or 



Eqs. 51, 67, 68 of Enfi lin & Voetl ( 120031) '). it has not been con- 
sidered that the mean RM is subtracted from the observed RM 
maps before our analysis is applied to them in order to remove 
contributions from foreground RM screens (e.g. the galaxy). 
We take this into account by noting that the RM value of a pixel 
at JCj. in the map will be changed only in the case of x± € Q., 
which can be written as 

RM(jc^) ^ RM*(;«:^) = RM(jc^)-l,.,en) f ^ RM(;«:1).(30) 

J An 

Using the properties of the Fourier transform of a function 
with compact support, one can show that 



where 



AnikJ 



(31) 



(32) 



and C^^{kj_) is the expected, uncorrected response to w(k). The 
linearity of the problem ensures that this correction is valid for 
any input power spectra w(k) (as long as C^^ik^) is calculated 
from it) which could be a delta-function, a power law, or could 
have any other shape. Similarly, we can write for the spheri- 
cally averaged expected observed power spectrum 



1 r^" 
— ld(/> 

^-^ Jo 



# |An(/t(cos 0, sin 0)r , (33) 



where w'^^'^ik) is given by Eq. ( I25I I. 

The subtraction of the mean RM does not only remove un- 
wanted homogeneous foregrounds, it also ensures that wT^{k) 
(or C'^^^,(kj^)) vanishes at the origin vvr''(0) = (similarly 



Cj^j^,(Oj_) = 0). For a correct behaviour of the responses at 
low k this correction is therefore crucial and will be taken into 
account in the following. 



3. Application to existing Rotation IVIeasure Maps 
of Galaxy Clusters 

Before we applied our analysis to the data any mean RM value 
was subtracted from the RM maps. This is necessary as men- 
tioned above, since we assume statistically isotropic fields and 
therefore we have to remove any RM foreground which would 
affect our analysis. Furthermore, very noisy pixels and pixel 
areas exhibiting huge differences in the RM values on scales 
smaller than the beamsize were removed from the data to de- 
crease the possible influences of observational artifacts. 

The origin of the z-axis was placed in the cluster centre and 
was used as reference point Xref for the global distribution of 
electron density and magnetic energy. The latter was assumed 
to scale with the electron density such that (B^(jr))'^^/Bo = 
{tieix) / rieo)"" ■ Unless stated otherwise, the parameter ob was 
chosen to be unity. For the electron density distribution a stan- 
dard yS-profile ^ was used. 

3.1. 3C75, 3C465 & Hydra A - the data 

The radio source 3C465 (or NGC 7720) is located in the 
Abell cluster 2634 of richness cla ss I. The redshift of the ob- 
ject is 0.0302 llWegner et al. 1999). A detailed X-ray study of 
this cluster was performed bv Schindler & Prieto (.1997) using 
ROSAT PSPC data. They derive for the density profile a core 
radius of 4.9 arcmin, ayS of 0.79 and an electron density at the 
cluster centre of «eo = 0.0012 cm"^. Their analysis revealed 
indications of a weak cooling flow in the cluster centre. The 
scenario of a weak cooling flow in Abell 2634 is also supported 
bv White (.2000) who analysed ASCA X-ray data. 

The radio galaxy 3C75 which is in the centre of Abell 
400 also of richness cla ss I has a redshift of 0.02315 
dde Vaucouleurs et alJll99l |). The X-ray properties of the clus- 
ter have been studied bv iRei prich & Bohrinser ( 200^) using 
ROSAT PSPC data. They determine for the gas density profile 
a core radius of 3.9 arcmin, a y6 of 0.54 and an electron den- 
sity at the cluster centre of Ugo = 0.0016 cm~^. There are no 
indications of a cooling flow in Abell 400 (White 2000). 

The details of the radio data reduction of the RM maps for 
the two sources above can be found in Eilek&Owen (20(jl. 
The typical RM values observed for them ranging between 
-250 rad m"^ and 250 rad m"^. The beamsize for the map of 
3C75 (Abell 400) is 3 arcsec which translates into 1 .4 kpc and 
for 3C465 (Abell 2634) the beamsize is 3.8 arcsec which is 
equivalent to 2.3 kpc for the cosmology chosen. 

Furthermore, it was assumed that the source plane is par- 
allel to the observer plane since no indication for the Laing- 
Garrington effect was found. 

Hydra A (or 3C218) is an extende d extragalactic radio 
sourc e located at a redshift of 0.0538 dde Vaucouleurs et alJ 
Il99ll) in the cluster Abell 780 of richness class 0. However, 
we will refer to it hereafter as Hydra A cl uster. Detailed X- 
ray studies have been performed on it (e.g. Ilkebe et al]ll997t 
IPeres et aljri998l bavid et alJl200ll) . For the derivation of the 



- defined as ng{r) - n^i^ + (r/rcf') where is called the core 
radius. 
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electron de nsity profile parameter, we relied on the work by 
lMohretal.l (fl999) done for ROSAT PSPC data while using the 
deprojection of X-ra y surfac e brightness profiles as described 
in the Appendix A of lPfrommer & Enfiliai2003 ). Since Hydra 
A is known to exhibit a strong cooling flow as observed in the 
X-ray studies, we assumed a double y6-profile ^ and used for the 
inner profile nei(O) - 0.056 cm"^ and = 0.53 arcmin, for 
the outer profile we used 11^2(0) - 0.0063 cm"-* and - 2.7 
arcmin and we applied a/3 - 0.77. 

The details of the radio data reduction can be found in 

IXav lor & Perlev ( 1993). The beamsize of the map is 0.3 arc- 
sec translating into 0.3 kpc. The source consists of two lobes, 
a northern and a southern one. Typical RM values in the north 
lobe are in the range between -1000 rad m"^ and -h3300 rad 
m"^ whereas in the south lobe also values of down to -12000 
rad m"^ were observed. 

We concentrated our analysis on the north lobe of this 
source because the signal-to-noise ratio in the data of the 
south lobe does not seem sufficient enough for our purposes. 
Furthermore, the application of the data filter to remove big 
steps between data pixels leads to a splitting of the RM data in 
the south lobe into a lot of small spatially disconnected areas. 
Such a fragmented window function can heavily obscure any 
power spectrum measured. 

For Hydra A, there is a clear depolarisation asymmetry of 
the two lobes observed as described by the Laing-Garrington 
effect. Taking this effec t into a ccount and following the anal- 
ysis of iTavlor & PerlevI ( 1 19931) . we first rotated the coordinate 
system about an angle of 1 10 degrees such that the x-axis was 
parallel to the jet axis. Afterwards the source plane was rotated 
about the new y-axis by an inclination angle 9 of 45 degree 
such that the north lobe would point towards the observer. 

3.2. Real Space Analysis 

After the mean RM subtraction and the noisy pixel filtering, 
the RM autocorrelation function Crm('"±) was determined us- 
ing Eq. (|5} choosing the described normalisation scheme. The 
resulting correlation function for Abell 2634, Abell 400 and 
Hydra A are exhibited in Fig. ^ where the correlation func- 
tion of Hydra A was divided by 100 for better display. The 
integration of these functions employing Eq. il2\ leads to the 
information on the RM correlation length Arm which was de- 
termined to be 7.9 kpc for Abell 2634, 5.3 kpc for Abell 400, 
and 1 .9 kpc for Hydra A. 

The determination of the magnetic autocorrelation function 
in real space w(r) achieved by the deprojection of the RM au- 
tocorrelation is numerically difficult due to the term in the de- 
nominator ^J)P■ - of Eq. and Eq. J10> which represents 
an integrable singularity and can in principle be avoided by a 
well selected coordinate transformation. 

However, the numerical calculation, especially the deter- 
mination of the value for r = of the magnetic autocorrelation 
function w{r) is difficult because an extrapolation to r = of 
the function w{r) being itself subject to extrapolation and data 

3 defined as n^r) = [nj,(0)(l + (r/r,,fr"' + /4(0)(1 + 




0.1 1 10 100 

r [kpc] 



Fig. 1. RM autocorrelation function |Crm('"±)I of Abell 400, 
Abell 2634 and Hydra A obtained using Eq.(|5}. Note that the 
correlation function for Hydra A was divided by 100 for repre- 
senting purposes. 



binning is involved. Therefore the determination of the point 
w(r - 0) is not accurate and one should be careful with the 
interpretation of the values derived by this method. The be- 
haviour of the magnetic autocorrelation around r = is exhib- 
ited in Fig.|2for the case of Abell 2634, where the magnetic 
field autocorrelation functions w{r) derived by employing vari- 
ous methods are shown. One can see that the steeper the initial 
slope of the function the less precise it becomes to determine 
w(r = 0) since slight deviations in the slope can lead to very 
different results. 

One should keep this in mind when calculating the value 
of the magnetic field autocorrelation length Ab which was de- 
termined to be 6.0 kpc for Abell 2634, 3.9 kpc for Abell 400 
and 1 .4 kpc for Hydra A using Eq. il4l . where for the value of 
w(r = 0) = By the extrapolation of Eq. to r = was used. 
Thus, the field strengths were determined to be of about 4 /jG 
for Abell 2634, of about 9 //G for AbeU 400 and of about 12 
juG for Hydra A in the real space analysis. 

3.3. Fourier Space Analysis 

For the calculations in Fourier space a 2-dimensional Fast 
Fourier Transform (FFT) algorithm was employed while set- 
ting all blanked values in the original RM map to zero. The RM 
autocorrelation function C±(k±) was then obtained by summing 
over rings in k-space. Since the magnetic field autocorrelation 
function w(k) is related to the RM autocorrelation function in 
Fourier space simply by multiplying the RM correlation func- 
tion C_i_{kj} by two (Eq. (I16» . there is no numerical difficulty 
involved deriving this function and thus, the results for the 
magnetic field quantities are more precise than those obtained 
applying the real space analysis. 

The integration of the power spectrum w°^^(k) following 
Eq.(ll9ll and Ea.( l20> yields the correlation lengths. Thus, the 
RM autocorrelation length Arm was calculated to be 8.0 kpc 
for Abell 2634, 5.3 kpc for AbeU 400 and 2.0 for Hydra A. 



8 



Corina Vogt and Torsten A. EnBlin; Measuring Cluster Magnetic Field Power Spectra 



20 

18 
16 
14 
12 
% 10 

f ' 
6 

4 

2 



-2 



using Eq.(IO) 
using Eq.(9) 
using w(l<) 



10 
r [l<pc] 



15 



20 



0.015 



0.01 



0.005 



-0.005 



-0.01 



-0.015 



"I I I I I I" 

using Eq.(IO) 

using Eq.(9) 
using w(k) 




..I I I.. 



20 40 60 80 100 120 140 160 180 200 
r [kpc] 



Fig. 2. The magnetic autocorrelation function w(r) for Abell 2634 calculated using Eq. (|9j, (I1Q> and i2l\ . One can clearly see 
that the value of w(r = 0) can deviate by a factor of two but for larger r all methods yield almost the same results. 



The magnetic field correlation length Ab was determined to be 
4.0 kpc for Abell 2634, 2.3 kpc for Abell 400 and 0.5 kpc for 
Hydra A in the Fourier analysis. 

No numerical problem is involved in the determination of 
the value of the magnetic autocorrelation function w(r) for r = 
if one uses Eq. <2H resulting in magnetic field strengths of 
about 5 yuG for Abell 2634, of about 1 1 fiG for Abell 400 and 
of about 23 juG for Hydra A. 

The application of Eq. J21> to the data translates the mag- 
netic autocorrelation function in Fourier space to real space. 
The comparison of the so derived function with the deprojected 
functions w{r) in real space is exhibited in Fig.|2lfor the case 
of Abell 2634. One can clearly see that at the origin there are 
deviations by a factor of two due to the uncertainties connected 
to the extrapolation used in the real space approach. At higher 
separations r all three functions shown do not differ signifi- 
cantly from each other. It is remarkable that except at the origin 
all three independent approaches are in such good agreement 
which demonstrates reliable numerics. 

The knowledge of the 3-dimensional power spectrum w{k) 
also allows calculating the magnetic energy spectrum s°g^{k) by 
employing equation J17> . The results for the three clusters are 
shown in Fig.|3l The magnetic energy spectra are suppressed at 
small k by the limited window size and the subtraction of the 
mean RM (small k in Fourier space translate into large r in real 
space). A response analysis as suggested above is performed in 
Sect. l3.4l in order to understand this influence in more detail. 

Another feature of these energy spectra is the increasing en- 
ergy density at the largest ^-vectors in Fourier space and thus, 
small r in real space. They can be explained by uncorrected 
noise on scales smaller th an the beamsize whi ch possesses a 
power spectrum of es oc ^ jEnBlin & Vogi2003l) . The determi- 
nation of the increasing slope of the energy densities for large 
k revealed proportionalities of eb ^' O- '-^ also shown in Fig. 
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Fig. 3. Magnetic energy spectrum s°i^^{k) derived for Abell 
400, Abell 2634 and Hydra A. The Fourier transformed beam- 
sizes (kbeam - ^Iheam) are represented by the thick vertical 
lines. The other straight lines describe the slope (dashed line 
represents sb °^ k and the dotted represents bb k k^'^) of the 
increase in the energy density for the largest k. 



13 Thus, we conclude that at the largest ^-scales uncorrected 
noise dominates the shape of the energy spectra. 

It seems reasonable to introduce an upper cutoff for the in- 
tegration of the magnetic energy spectra in Eq. Ml\ . In Fig. 
|3j the equivalent beamsize in Fourier space kbeam - ^/Ibeam 
(where Iheam is the beamsize in real space defined as FWHM) is 
represented by a vertical line for each cluster which is 1 .4 kpc"' 
for Abefl 2634, 2.2 kpc ' for Abell 400 and 10.0 kpc ' for 
Hydra A. One can clearly see that the noise induced increase 
of the energy density hes on ^-scales beyond kbeam- Therefore, a 
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suitable upper cutoff for any integration of the magnetic energy 
spectrum seems to be kheam- 

However, the influence of the upper A:-cutoff in the integra- 
tion can be seen in Fig. which displays the magnetic field 



10 



strength B = y/^nisBj estimated from 

r 

Jo 



SB(k < kc) 



dkss(k). 



(34) 



100 




100 



Fig. 4. Magnetic field strength B dependent on the upper k- 
cutoff in the integration of Eq. jilt . The equivalent beamsizes 
kheam are represented as vertical lines. 



3.4. Test the model and assess the influence of the 
window function 

The various possibilities to test the window function and to as- 
sess its influences on the results outlined in section IZTI were 
also applied to the data. The ;i'^-function was derived by radi- 
ally averaging the x^(^±) distribution resulting from Ea. (l22> . 
where for the magnetic autocorrelation length Ab the value 
derived in the Fourier analysis was used since this approach 
seems to be more accurate. The resulting radially binned dis- 
tributions are shown in Fig.|5]determined for the three clusters. 
There is no apparent spatial large scale trend visible for Abell 
400 and Abell 2634 which indicates a reasonable model for the 
window function. In the case of Hydra A, there appears to be a 
trend of higher values for x^(x±) towards larger r, comparable 
to small scale trends seen in the x^i^^) distribution of Abell 
2634 and Abefl 400. 

Furthermore, thex^i^±) distribution was integrated follow- 
ing Ea.( l23t and;^'^,, was calculated to be 1 .0 for Abell 2634, 1 .2 
for Abell 400 and 1 .6 for Hydra A. A refinement of the model 
describing parameter does not seem to be required for Abell 
2634 and AbeU 400. 

Since the value forxl^ of 1.6 for Hydra A is deviating from 
1, we varied the geometry of the source, i.e. we varied the in- 
clination angle 9 between source and observer plane from 30 
degree up to 60 degree resulting in field strength ranging from 
17 nG to 33 juG and values for xiv of 1-4 to 1.8, respectively. 
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Fig. 5. Testing the window function by calculating the x^(x±) 
distribution for all three cluster 



These values were derived by integrating over the full accessi- 
ble unrestricted (without cutoffs) A:-space. 

Furthermore, we varied the scaling parameter ag in the 
scaling relation of the electron density to the magnetic energy 
density. A value for ub of 0.5 is still in the li mit of reason- 
able values as suggested by bolaget alJteOOll) . Furthermore, 
such a scaling parameter of aB = 0.5 means that the magnetic 
energy sb{x) oc ndx) and thus, the magnetic field would be 
proportional to the thermal energy of the cluster gas assum- 
ing approximately isothermal conditions. However, in this case 
(ob = 0.5) one obtains for xl^, a value of 1.2 and the magnetic 
field strength is reduced to 17 fiG by integrating over the full 
accessible ^-space. 

One can not be sure that the cause of the trend in thex^ix±) 
distribution is due to a geometry other than assumed because it 
could also be explained by a fluctuation similar to the one seen 
in the^^(xj.) distributions of Abell 2634 and Abell 400. There 
is no reason to change the initial assumption for the geometry 
of Hydra A. However, we should keep in mind that the central 
field strength given could be slightly overestimated for Hydra 
A but it is not clear to what extent if at all. 

A good test for the validity of the isotropy assumption is 
the inspection of the Fourier transformed RM data |RM(fc)p as 
shown in Fig.|6lfor Abell 400 and Hydra A. The EFT data look 
similar for Abell 2634. In this figure half of the Fourier plane 
is shown since the other half is symmetric to the one exhibited. 
No apparent anisotropy is present in this figure. 

We mentioned already that the magnetic energy spectrum 
s°^'^(k) in Fig.|3]is suppressed by the limited window size for 
small fe-vectors in Fourier space. Therefore a response analysis 
is necessary in order to understand the influence of the window 
on the shape of the magnetic energy spectrum. For this pur- 
pose, the projected window function was Fourier transformed 
by employing a EFT algorithm and inserted into Eq. i25i . We 
then compared the response functions vv'^''P(A:) obtained to the 
observed power spectrum w°^'^(k) as shown in Fig. 0for the 
case of Abell 2634, where we chose for the normalisation a 
magnetic field strength Z? of 5 yuG. The corresponding figures 
look very similar for Abell 400 and Hydra A. 
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Fig. 6. The power \(SM(k)}\^ of the Fourier transformed RM 
map of Abell 400 on the left side and Hydra A on the right 
side. The image for Abell 2634 looks similar. Only half of the 
Fourier plane is exhibited the other half is point symmetric to 
the one shown 
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Fig. 7. Responses at various scales p to the window function of 
Abell 2634 assuming a B = 5//G in comparison to the observed 
magnetic autocorrelation function w(k). 

From Fig.Q one can clearly see that the response to delta 
function like input power spectra on small p scales in Fourier 
space (i.e. large r in real space) is a smeared out function as one 
would expect. The response for larger p becomes more strongly 
peaked suggesting that at A;-scales larger than 0.3 kpc"' for 
Abell 2634 the influence of the window function becomes neg- 
ligible. From similar plots for the other two clusters under con- 
sideration, scales of about 0.4 kpc"' for Abell 400 and 1.0 
kpc"' for Hydra A are found. Thus, one would use this value as 
a lower cutoff k^^i,^ for the determination of the magnetic field 
strength. 

From the response power spectra Wp{kj_) calculated for dif- 
ferent p requiring a magnetic field strength B, one can derive a 
magnetic field strength Bexp by integration. In Fig.|8j the com- 



parison between these two field strengths is plotted for the dif- 
ferent /7-scales. It can be seen that for the smallest p's the devi- 
ation between expected field strength and actual observed one 
is significantly but for larger p they are almost equal suggest- 
ing that on these scales the influence of the window is negligi- 
ble. Taking Fig.|Slinto account, the values determined above as 
lower cutoff k^in can be confirmed. 
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Fig. 8. Comparison between the expected magnetic field 
strength B to the measured one Bgxp for the particular p-scales 
for the three clusters. 



Another possibility to assess the influence of the window 
on the power spectra is to match the added response functions 
and the actual observed power spectra w°^^(k). This was done 
by applying Eq. ( I27t to a set of response functions generated 
for closely spaced /^-scales. The resulting function w^^'^'ik) was 
matched to the actual observed power spectrum w°^^{k) for each 
of the three clusters by varying Bq, the spectral index a of the 
power law in Eq. J27> to match the slope of the functions and 
Pmin as lower cutoff to fit the function at its turnover for small 
^-scales. The resulting functions are shown in Fig.|9]in compar- 
ison to the respective observed power spectra. The shape of the 
power spectra was matched for Abell 2634 for a spectral index 
of a = 1.6, for AbeU 400 for an a = 1.8 and for Hydra A for 
an ff = 2.0 whereas a lower ^-cutoff of 0.09 kpc"' was used for 
Abell 2634, 0.08 kpc"' for Abell 400 and 0.3 kpc"' for Hydra 
A. 

The spectral index a can also be used to plot the respec- 
tive power law as represented by the straight lines in Fig. |9] 
The slopes of these lines clearly deviate from the respective ob- 
served power spectra suggesting that one cannot estimate dif- 
ferential parameters like spectral indeces directly from Fourier 
transformed RM maps. More sophisticated approaches should 
be developed. However, our approach allows to exclude flat- 
ter spectral slopes than a - 1.3 still leaving a Kolmogorov 
spectrum (a = 5/3) as a possible description. For attempts to 
measure the magnetic power spectrum from cluster simulations 
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Fig. 9. Closely /7-spaced response functions were integrated 
by employing Eq. and the resulting M>'^''P(fc)are exhibited in 
comparison to the observed power spectrum ■w°'°^(K) shown as 
solid lines. The spectral index a used for Abell 2634 is 1.6, for 
Abell 400 is 1.8 and for Hydra A is 2.0. The straight lines rep- 
resent the power law assuming the spectral index determined 
from the response analysis. Note that the data for Abell 2634 
are multiplied by 0.1 and the data for Hydra by 10 for repre- 
senting purposes. 

and radio maps see lOolae et alJ ( l2002l) and lGovoni et alJ (l2002h 

respectively"*. 

4. Results 

A summary of the values for RM autocorrelation length /Irmi 
magnetic field autocorrelation length and the central mag- 
netic field strength Bo for the magnetised cluster gas in Abell 
2634, Abell 400 and Hydra A, derived in Fourier space and real 
space is given in Tab. ^ The RM autocorrelation length /Irm 
determined in both spaces is almost equal which demonstrates 
reliable numerics. Furthermore, it is found that the RM auto- 
correlation length /}rm is larger than the magnetic autocorrela- 
tion length Xb- The deviation of the magnetic field quantities 
between both spaces is caused by the numerically complicated 
deprojection of the magnetic autocorrelation function w(r) in 
real space. Thus, the results in Fourier space are more reliable 
and should be considered for any interpretation and discussion. 

In Sect. 13.41 we discussed many possible influences on the 
power spectra. We verified that for intermediate A;-ranges ex- 
tending over at least one order of magnitude the power spectra 
are not governed by the window or the resolution of the RM 
maps and thus, represent most likely the magnetic field prop- 
erties of the cluster gas. Thus, one can derive field quantities 
such as magnetic field strength Bq and autocorrelation length 
Xb from this intermediate A:-range. 

The suitable cutoffs for the integration of the power spec- 
tra in order to derive the magnetic field strength and autocor- 
relation length were determined in Sec. 13.41 Using for Abell 

The conventions describing the spectra may differ in these articles 
from the one used here. 



2634 a lower cutoff fc,„,„ of 0.3 kpc"' and an upper cutoff of 1.4 
kpc ' results in a magnetic field strength of 3 /vG. Integrating 
the power spectrum of Abell 400 in the range from 0.4 kpc"' 
to 2.2 kpc"' yields a field strength of 6 /iG. The same was done 
for the Hydra A cluster in the range from 1 kpc"' to 10 kpc"' 
which resulted in a field strength of 13 yd. The determination 
of the magnetic field autocorrelation length Xb by integration 
over the power spectra in limited A;-space yields 4.9 kpc for 
Abell 2634, 3.6 kpc for Abell 400 and 0.9 kpc for Hydra A. 

In the course of the response analysis in order to esti- 
mate the influence of the window on the power spectra, we 
matched the integrated response power spectra with the 

actual 3-dimensional observed power spectra vv°''^(A;) and de- 
termined the spectral index o-, cq and the lower cutoff Pmin in 
Eq. (|^). Using these values for the parameters in Eq. J28> , 
the direct integration in the limits of and /?,nax = ^beam of 
the power spectra obtained by the response analysis results in 
magnetic field strength of about 3 /iG for Abell 2634, of about 
6 iiG for AbeU 400 and of about 1 1 //G for Hydra A. These 
field strengths are in good agreement with the results obtained 
by the restricted Fourier space integration of the observed 3- 
dimensional power spectra. 

However, the values for Xb were determined for the same 
Pmin and pmax = ^beam applying Eq. (|29|l. The resulting value 
is 13 kpc for Abell 2634, 17 kpc for Abell 400 and 4 kpc for 
Hydra A. They deviate by a factor of 4 to 5 from the mag- 
netic field autocorrelation length Xb determined from the re- 
stricted Fourier space integration. The reason for this behaviour 
could be found in the sensitivity to the lower cutoff /?min and 
the power law index a. We have already mentioned that our 
method might not be suitable at that stage to determine differ- 
ential parameters such as a power law index a. A more detailed 
discussion which is beyond the scope of this paper is necessary 
in order to understand this behaviour. 

However, the values of the magnetic field strength Z?o for 
Abell 400 and Abell 2634 of about 6 yuG and 3 yuG, respec- 
tively, are larger comp ared to the initial analysis performed by 
lEilek & OwenI ( l2002l) for the case of a magnetic sheet with a 
thickness of 10-20 kpc. They discuss a variety of magnetic field 
models and we refer to Tab. 2 of their paper for detailed num- 
bers and description. 

For the Hydra A cluster, iTavlor & P erlev' ('l993*) conclude 
in their analysis that there is a random magnetic field compo- 
nent of about 30 yuG with a correlation length of 4 kpc. The 
deviation by a factor of 3 of the field strength from the value 
revealed by our analysis might be explained by the usage of 
an improved electron density profile for our analysis which 
also takes the cooling flow into account. Our value for the 
central magnetic field strength Bo might be also lower due to 
the conservative approach of restricting the ^-space integration 
range. Another explanation for the difference could be sought 
in our exclusion of the south lobe from the calculation above. 
Including the south lobe in our analysis leads to higher central 
field strength but given the influence of the very complicated 
window function in the case of the south lobe, it is not clear to 
what extent the real power spectrum is resembled. 

The estimation of the dynamical importance of the mag- 
netic fields derived for the cluster gas can be done by compar- 
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Real Space 




non-restricted Fourier Space 


restricted Fourier Space 


cluster 


cs 


/Irm [kpc] 


Ag [kpc] 


Bo LuG] 


/Irm [kpc] 


Ab [kpc] 


Bo UiG] 


Ab [kpc] 


Bo ImG] 


Abell 400 


1.0 


5.2 


3.9 


8 


5.3 


2.3 


12 


3.6 


6 


Abell 2634 


1.0 


7.9 


6.0 


4 


8.0 


4.0 


5 


4.9 


3 


Hydra A 


1.0 


1.9 


1.4 


13 


2.0 


0.5 


23 


0.9 


13 


Abell 400 


0.5 


5.2 


3.9 


7 


5.3 


2.3 


10 


3.4 


6 


Abell 2634 


0.5 


7.9 


6.0 


3 


8.0 


4.0 


4 


4.9 


3 


Hydra A 


0.5 


1.9 


1.4 


10 


2.0 


0.5 


17 


0.9 


9 



Table 1. Values for the autocorrelation length scales Ab and Arm and the magnetic field strength Bq at the cluster centre obtained 
in real space analysis, in Fourier space analysis and in restricted Fourier Space analysis are given for the different clusters under 
consideration. The values are calculated for two different scaling parameters a^. Note that for position other than the cluster 
centre the average magnetic energy is given by {B^{x)y^^ - Bo (neix)/neo)"'" . 



ing the thermal pressure (p,h - 2ne{Q)kT core) with the mag- 
netic pressure {pB = B^H^n)). One can calculate pbI p,h which 
yields 0.08 for the case of Abell 2634 (assuming a Tore -1-2 
keV (Schindler & Prieto 1997)), pb/pu, = 0.19 for Abell 400 
(Tcore = 1-5 keV (lEilek & O wen 2002)) and pr/p,,, = 0.01 for 
Hydra A (Tcore ^ 2.7 keV ( jPavid et aLi200ll) '). It is astonish- 
ing that the value of PbIpih is smaller for Hydra A, which is a 
cooling flow cluster, than for the non-cooling flow cluster Abell 
400 and Abell 2634. The values of pbIpih = 0.1. ..0.2 for the 
latter two clusters give an indication that for those clusters the 
magnetic field is of some weak dynamical importance for the 
cluster gas. 

5. Conclusions 

We presented the application of a new analysis of Farad ay ro- 
tation maps recently developed bv lEnBlin & Vogll(l2003b to ob- 
servational data in order to estimate magnetic field strength and 
autocoiTelation length and in order to determine the magnetic 
power spectra of the magnetised intra-cluster gas. We described 
that the analysis relies on the assumption that the magnetic 
fields are statistically isotropically distributed throughout the 
Faraday screen. We introduced the window function through 
which any virtually statistically homogeneous magnetic field 
can be thought to be observed. This window function describes 
the geometry of the source and the global properties of the 
intra-cluster gas such as the electron density, known from X- 
ray measurements, and the global average magnetic field dis- 
tribution, which we assume to scale with the electron density. 
Furthermore, we explained two possible approaches in real and 
Fourier space and outlined the tests for the evaluation of any in- 
fluence especially arising from the observational nature of the 
data such as limited source size, resolution and pixel noise on 
the results obtained. However, we stated that our analysis al- 
lows to measure average magnetic energies but it is not sensi- 
tive to the particular realisation of the magnetic field structure. 

We applied this approach for the first time to observational 
data and derived not only reliable results for magnetic field 
strength and autocoiTelation length but also we got insight into 
strength and quality of any results obtained by applying this 
new analysis to any observational data. In order to understand 
the possible impacts on the results we reanalysed the Faraday 
rotation maps of three extragalactic extended radio sources, i.e. 



3C75 in the cluster Abefl 400, 3C465 in Abefl 2634, which 
were kindly provided by Frazer Owen and Jean Eilek, and 
Hydra A in Abell 780, which was kindly provided by Greg 
Taylor. We performed the analysis in real and Fourier space for 
these three Faraday rotation maps. While discussing the diffi- 
culties involved in the application to the data, we realised that 
the calculations in Fourier space are more reliable. 

We tested the isotropy assumption and no indication of 
anisotropy was found. Furthermore a ;t'^-test was performed in 
order to assess the model adopted for the geometry of the sam- 
pling volume incorporating the global electron density and av- 
erage magnetic energy distribution. In the case of Abell 2634 
and Abell 400, no indication was found that our model may 
be incoiTect. However, in the case of Hydra A we found in- 
dications that the window function needs refinement but the 
indications were not strong enough in order to enforce this re- 
finement. 

We realised that the magnetic energy spectra EB{k) of the 
three clusters investigated are dominated on the largest k- 
scales (i.e. smallest r-scales) by uncorrected noise. Therefore 
it seemed natural to introduce a higher A;-cutoff fc^ax for any 
integration in A:-space necessary to derive magnetic field quan- 
tities such as field strength and autocorrelation length. Being 
conservative, we used as higher A;-cutoff the equivalent beam- 
size in Fourier space k^^^, - n/lhgam which is equal to 1.4 kpc"' 
for AbeU 2634, 2.2 kpc ' for Abell 400 and 10 kpc ' for Hydra 
A. This will have the effect of loosing some power which is 
redistributed due to the window from smaller to these larger 
fc-scales. 

On the smallest A:-scales, i.e. the largest r-scales, power is 
suppressed because of the limited size of the window. For the 
assessment of the influence of this suppression on the power 
spectra, we applied a response analysis described in Sec. 12.41 
to the observational windows. The response of the window 
to delta like input power spectra on small p-scales is a wide, 
smeared out function whereas the response on larger p-scales 
is a peaked function. This observation motivated the introduc- 
tion of a lower A:-cutoff k^in in any integration in A:-space. This 
value was determined to be 0.3 kpc"' for Abell 2634, 0.4 kpc"' 
for Abell 400 and 1.0 kpc"' for Hydra A. However, the mag- 
netic field strengths derived are not sensitive to this lower cutoff 
due to the small k power suppression. 
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We verified that the intermediate A:-ranges between fc,ni„ and 
^max extending over at least one order of magnitude can be used 
to determine actual magnetic field properties of the intra-cluster 
gas. 

Integrating over the response power spectra on particular 
;?-scales enables us to match the so determined power spectra 
and the actually observed power spectra in order to have an 
independent estimate for magnetic field strengths and spectral 
slopes. For the three clusters under consideration, spectral in- 
dices a for the slopes of the power spectra were determined to 
be in the range from 1 .6 to 2.0. Therefore it would be possi- 
ble that the magnetic field in these clusters have a Kolmogorov 
power spectrum exhibiting a - 5/3. However, presently we can 
not exclude steeper spectra but flatter spectra exhibiting slopes 
smaller than or = 1 .3 would have been recognised by our anal- 
ysis. Although we realised that this analysis is not suitable for 
flie determination of differential parameters such as the spec- 
tral slopes of power spectra directly from Fourier transformed 
RM maps, the determination of integrated quantities such as 
the magnetic field strength appears feasible. 

Taking all these arguments into account we determined by 
integration of the magnetic energy spectrum in the limits be- 
tween kjnin and ^max values for the RM autocorrelation length 
/Irm, the magnetic field autocorrelation length Ab and the cen- 
tral average magnetic field strength Bq. An overview of these 
numbers is given in Tab. where these values are compared to 
those derived in real space and non-restricted Fourier space. 

The magnetic field autocorrelation length Ab was deter- 
mined for the restricted Fourier space integration to be 4.9 kpc 
for Abell 2634, 3.6 kpc for AbeU 400 and 0.9 kpc for Hydra A. 
In comparison with the RM autocorrelation length /(rm calcu- 
lated to be 8.0 kpc for Abell 2634, 5.3 kpc for Abell 400 and 
2.0 kpc for Hydra A, it can be said that these two characteristic 
length scales differ from those often assumed in previous work. 

The magnetic field strength in the cluster centre Bq was cal- 
culated for the same limited fc-space and was determined to be 
3 juG for Abell 2634, 6 /jG for Abell 400 and 13 fiG for Hydra 
A. Given the assumption of isotropy and a scaling parameter 
Ob - 1.0, these are conservative values. The resulting mag- 
netic pressures suggest some small but non negligible pressure 
support for a dynamical influence in the case of Abell 2634 
and Abell 400 since the value of PbIpiIi is 0.08 and 0.19, re- 
spectively. 

We note, that our approach so far is not able to separate 
the influence of noise on the maps from the astrophysically in- 
teresting signal of intra-cluster magnetic fields. Investigations 
aiming to d etect and minimise such spurious signals are in 
preparation llEnBlin et al ] il2003l) : Dolag et al., in prep.). 
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